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The density matrix renormalization group is applied to a relativistic complex scalar field at finite 
chemical potential. The two-point function and various bulk quantities are studied. It is seen 
that bulk quantities do not change with the chemical potential until it is larger than the minimum 
excitation energy. The technical limitations of the density matrix renormalization group for treating 
bosons in relativistic field theories are discussed. Applications to other relativistic models and to 
nontopological solitons are also suggested. 
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I. INTRODUCTION 

In relativistic field theory on the Euclidean lattice, the 
chemical potential manifests itself as an imaginary vec- 
tor potential in the imaginary time direction |l| . The use 
of Monte Carlo techniques and reliance on importance 
sampling has the effect of making the probability weight 
meaningless if the action or the fermion determinant is 
complex; this is known as the "sign problem" . With lat- 
tice fermions, the sign problem is also a consequence of 
a complex fermion determinant resulting from the inclu- 
sion of a chemical potential term. A variety of techniques 
are commonly used to circumvent this problem, includ- 
ing Glasgow reweighting Q and analytic continuation to 
an imaginary chemical potential Q. 

The aim of this paper is to use the density matrix 
renormalization group (DMRG) to study a toy model in 
relativistic quantum field theory with a nonzero chemical 
potential. The results of this paper demonstrate that 
this approach correctly captures the phenomenology of 
the theory. We seek the ground state of the Hamiltonian 
for the system using variational methods and so DMRG 
does not suffer from the aforementioned sign problem Q ■ 
This problem of relativistic field theories at finite density 
should be contrasted with the negative sign problem of 
fermions in quantum Monte Carlo [5|. It is already widely 
acknowledged that this is avoided by DMRG [g, 0| • As 
far as we are aware, however, no attempt has been made 
to circumvent the relativistic finite-density sign problem 
with DMRG. 

When the chemical potential \x is smaller than the low- 
est excitation energy (at zero coupling, the bare mass m), 
the bulk quantities studied are seen to have the same 
value as at fi — 0. This effect is known from QCD as 
the "Silver Blaze" problem [8|. Such an effect cannot be 
studied in a nonrelativistic field theory. 

For larger chemical potentials, the system condenses. 
We cannot voyage far into the condensed phase because 
the truncation in states converges more slowly. For 



smaller values of the chemical potential, however, conver- 
gence is better. Nonetheless, the truncation in bosonic 
states and the limitation to (1 + 1) dimensions remain 
the main drawbacks of using DMRG as a nonperturba- 
tive tool in quantum field theory. 

The density matrix renormalization group has previ- 
ously been applied to bosonic problems with nonzero 
chemical potential for condensed- matter systems [9|. In 
such cases the models usually describe a single bosonic 
field and there is no symmetry between the particles and 
antiparticles. Such Hamiltonians emerge as the nonrel- 
ativistic limit of the model considered here when the 
parts of the Hamiltonian concernin g th e antiparticles 
and their interactions are neglected The current 

work, however, accurately describes the fully relativis- 
tic U(l) model. Toy models in particle physics have 
been studied using DMRG previously, in particular a one- 
component scalar field with A</> 4 interaction [llj , the mas- 
sive Schwinger model 1121 and a simple model exhibiting 
asymptotic freedom [13j . SU(iV) spin chains have also 
been studied [l4j |. 

In this paper we consider a scalar field model with 
U(l) symmetry in (1 + 1) dimensions. For convenience, 
we transform fields to a two-component real scalar field, 
giving a Lagrangian 
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) - -m 4>n4>n ~ -\(4>n4>n) • (1) 



The conjugate momenta are n n — 4> n . We then perform 
a Legendre transform to give the corresponding Hamil- 
tonian density for this theory [15| . 



Ho = i(V0„)(V^„) 
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2 2 
This model has one conserved charge 



Q = J dx j = J dx [0l7T 2 - </> 2 7Tl] . 



(2) 



(3) 
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We introduce the chemical potential fi as a Lagrange mul- 
tiplier into an effective Hamiltonian rt for minimizing the 
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energy at nonzero total charge Q and obtain 
H=H - Mo 

= ^(W„)(Wn) + 7j7r„7T n + ^ m 2 <j) n (j) n + ~A(0 n ^ n ) 2 
- H (017T2 - 4>2^l) , 

(4) 

which we will use for our studies of this model using the 
DMRG. 

This paper is organized as follows. In Sec. |TT]we out- 
line the density matrix renormalization group as applied 
to the present model and discuss places where our im- 
plementation differs from those in the literature. Our 
numerical results are presented in Sec. IIII1 and our final 
remarks can be found in Sec. IIV1 



II. THE DENSITY MATRIX 
RENORMALIZATION GROUP 

The density matrix renormalization group is a varia- 
tional technique for finding quantum states of quasi-one 
dimensional systems (l6| . Originally conceived to study 
systems - such as the Heisenberg model - on lattices 
too large to treat with exact diagonalization, it is essen- 
tially a development of Wilson's numerical renormaliza- 
tion group for handling interacting systems ((J [l?], EH ■ 
A brief qualitative summary is given here. The entire 
system is termed the "superblock," and is divided into 
two renormalized "system blocks" , between which one or 
more sites are "inserted" . This allows numerical diag- 
onalization of a smaller Hamiltonian than if the entire 
system were exactly diagonalized. Typically, one is only 
interested in the ground state, and so only one eigen- 
vector needs to be numerically obtained. Operators are 
transformed to best represent the states of interest and 
the inserted site is incorporated into one of the system 
blocks. The process repeats with a new "inserted site" . 

When this technique is used to study a finite system 
as in the present case, the representation of the system 
is optimized by repeated "sweeping" to one end of the 
lattice and then the other. The system size is kept con- 
stant by growing one system block at the expense of the 
other. Sweeping stops when observables of interest no 
longer change. If the system is homogeneous then the 
symmetries of the system can be exploited to accelerate 
this process. 

Our approach is to insert only a single site, rather than 
two sites, during the sweeping process. This approach 
has been used successfully to study the one-component 
A0 4 model in Ref. fllj ]. and its use remains appropriate 
here. Generally, when fermionic systems are considered, 
the number of states per lattice site required to fully de- 
scribe the system is not large and indeed for small lattices 
it would in principle be possible to exactly diagonalize the 
system. In addition, there are unambiguous reasons why 
adding two sites helps to improve efficiency when using 



the infinite- volume algorithm with spin systems [l6j , in- 
cluding that it keeps the system symmetric. However, in 
our present system we must truncate the tower of bosonic 
states and even for a small system there is no way to ex- 
actly diagonalize the Hamiltonian. As each inserted site 
gets incorporated into our renormalized blocks the high- 
est state accessible is given by the highest state in the 
truncation used for the single site. 

Let us label the number of states kept in the system 
block by M and the number of states on the inserted 
site N. By the above argument it seems reasonable to 
expect that, for bosonic systems with a given maximum 
dimension D of the Hamiltonian, the best numerical re- 
sults will be obtained by taking the largest possible trun- 
cation for three sites (D ss M 2 N) rather than for four 
(D « Af 2 iV 2 ). 

Unlike in Monte Carlo simulations, where the prefer- 
ence is for periodic or twisted boundary conditions, for 
DMRG calculations the results are more precise when ob- 
tained with open boundary conditions. The interaction 
between the extreme ends of the lattice cannot be renor- 
malized within the traditional scheme for finite system 
DMRG and hence the system only has an approximate 
translational invariance. A translationally invariant ap- 
proach is possible when one reformulates the problem in 
terms of matrix-product states [ill [2(| ■ Treating bosons 
will remain difficult, due to the high bond dimension of 
the matrices needed. For the purposes of this paper, it 
is sufficient to consider the 'traditional' choice of boson 
number states, with open boundary conditions. 

A. The model 

A single-component relativistic scalar field was first 
studied with the DMRG in Refs. [U El- This discus- 
sion therefore parallels these works, extended to a two- 
component field. Moving to a (L x oo) lattice, we keep 
the conjugate momenta as operators 7r„ but discretize the 
gradient term in in the usual way, 

V<fi n ^-[(f> n {x + a)-(f> n (x)], (5) 
a 

where we have introduced a lattice spacing a and label 
the sites of the lattice by 0, . . . , x, x+a, . . . , L. We can set 
a = 1 without loss of generality as we can vary m, A and 
H instead. We work with real-space DMRG (although 
momentum-space formulations exist), which means we 
need a set of basis states for every lattice site. We treat 
4>\ and 4>2 as separate fields with their own creation and 
annihilation operators, 

<f>n(x) = (al(x) + a n (x)) , (6) 

n = 1,2, and 
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where a^ n and a n create and annihilate particles of type 
n at sites labeled by x on the lattice. This motivates the 
use of \p, q) for boson number states at a given lattice 
site, where p and q label each of the two different particle 
types. 

The equal time commutation relation, 

[(f)„(x),Tr m (y)] = i5 x , v S m ,m (8) 

becomes 

[a n (x),o^(3/)] = 5 Xj y5 m>n . (9) 
The Hamiltonian becomes 

+ \\[M*) 2 + M*f\ 2 ) ( 10 ) 
+ J2l(lMx)-Mx + i)} 2 

x=l 

+ [0 2 (x)-^ 2 (x + l)] 2 ) . 
We then arrive at the effective Hamiltonian on the lattice 

L 

H = H -J2v [iT2{x)Mx) - *i{x)fo{x)] . (11) 

x=l 

Quantized as outlined above, Hq can always be written 
as a real symmetric matrix. A nonzero chemical potential 
requires us to diagonalize a complex, Hermitian Hamil- 
tonian H but this has been done previously in DMRG 
studies of electron systems with persistent currents [22j . 
It does not present any problem for the DMRG. 

We split the Hamiltonian up into the Hamiltonians for 
the two renormalized 'system' blocks and the single in- 
serted site, plus interaction terms [ll| . 

H = Hl + h n -i in + h n + h„ tn+ i + Hr (12) 

where Hl and Hn are the left and right system blocks 
respectively. The single site h n is 

h n = ~ [Ki{n) 2 + TT 2 {n) 2 ] + i?7i 2 [(f)i{n) 2 + M™) 2 ] 

1 2 

+ -A [cf)i{n) 2 + 4>2{n) 2 ] +n [TT 2 (n)(j)i(n) - 7ri(n)<^ 2 (»)] > 

(13) 

and the interaction term ft. n , n +i is given by 

hn.n+l = g ([0lW-^l(«+l)] 2 + [02H-0 2 (n+l)] 2 ) ■ 

(14) 



Diagonalizing H numerically, we obtain an approxima- 
tion to the ground state \ip), 

\Tl>) = Y,i>nk\i)\j)\k), (is) 

where the single labels i, j and k run over the truncated 
bases. The left and right blocks are in general renor- 
malized (except at the ends of the lattice) and the basis 
states are optimized, but the central site's basis \j) al- 
ways corresponds to a sum over the two-particle basis 
states \p, q) . 

The matrix elements are stated for X(f> 4 with a single 
scalar field in Ref. (2lj . and have been summarized in a 
form generalized to the present case in Appendix El 

Regardless of what numerical technique we use to 
study the discretized Hamiltonian system, we will have 
to truncate the basis states in (|T5|) . This can be inter- 
preted as a UV cutoff that may not always be as high as 
that of the lattice spacing, I /a. 



B. Relation with microscopic Hamiltonians for 
Bose-Einstein condensation 

The density matrix renormalization group has been 
used in condensed-matter systems to study microscopic 
Hamiltonians that exhibit Bose-Einstein condensation, 
such as Bose-Hubbard models 0, [23| . These models are 
substantially different from the present case as there is 
no possibility of spontaneous symmetry breaking; there 
is only a single particle species present. In this section 
we show how the present work corresponds to relativis- 
tic physics that cannot be obtained with a microscopic 
Hamiltonian motivated by problems in condensed mat- 
ter. 

For the current discussion, let us work with free fields 
(A = 0). In condensed-matter systems, Bose-Einstein 
condensates are often modeled by the Gross-Pitaevskii 
equation. Following Ref. [Io|, we therefore want to trans- 
form to nonrelativistic fields "J and $ that (together with 
their complex conjugates) satisfy the Gross-Pitaevskii 
equations of motion with a nonrelativistic chemical po- 
tential fj, nr , 

(-i| + ^- /te )* = 0, (16) 

and equivalently for In terms of <pi and <f>2 (and their 
conjugate momenta), the appropriate transformation is 

Mi) + * = -^= (*(*) + **(*)), (17) 

7n(x) + i n 2 (x) = i^(**(x) - *(a:)) , (18) 

where ui = \Zk 2 + m 2 . One should then expand to order 
k 2 . If we compare the free-field action for the relativistic 
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fields 



2 

Srel = I dt I dx 

"* n=l 



d 2 d 2 

dt 2 dx 2 



+ m 2 )4> n , (19) 



with that for one of the nonrelativistic fields 



Sri 



dt / dx <r 



■ 9 1 „„ 



* (20) 



then we see that the nonrelativistic fields split into two 
parts, one for the $ field and one for the "J field. We have 
omitted the nonrelativistic fields ^ and ^* correspond- 
ing to one of the particle species, as they are suppressed 
by a factor e~^ +m '^ * in the partition function. The 
transformations (jl7H18[) have had the added effect of di- 
agonalizing the chemical potential term for each particle 
species - the term in the Hamiltonian referring to the 
chemical potential is now just /z 4 r *^ r , a number opera- 
tor that is usually combined with m to give /z nr . Hence, 
the sign problem disappears for nonrelativistic models at 
finite density as one can treat the system with one field. 



C. Extension to finite temperature 

Ideally, we would be able to study nonzero density 
at finite temperature. Since we keep the conjugate mo- 
menta 7r„ as operators throughout we are dealing with a 
(1 + l)-dimensional system, and cannot introduce finite- 
temperature physics in the standard, Euclidean, way. For 
the DMRG, the traditional method used to treat small 
T > is to find several low-lying states [f| [24| 



(21) 



when numerically diagonalizing the Hamiltonian, and 
weight them by the Boltzmann factor in the reduced den- 
sity matrix p to give 



P : 



-0E, 



E 



, (n) . * (n) 
Yijk'Vi'j'k'- 



(22) 



This adds an extra layer of truncation: we must make 
sure that the number of excited states included in the 
density matrix is large enough to correctly capture the 
finite-temperature behavior to the accuracy allowed by 
our truncation. The most straightforward way to imple- 
ment this is to vary the number of eigenvectors obtained 
by our numerical diagonalization code for inclusion in 
Eq. [5H If increasing the number of eigenvectors ob- 
tained does not affect the observables of interest then 
we can consider the truncation adequate. This ought to 
happen when the smallest obtained Boltzmann weight, 
associated with the highest energy, e~ l3Erl ^ x is compa- 
rable in size to our algorithm's convergence tolerance. 

These added complications mean that finite- 
temperature behavior for a bosonic system cannot 
be reliably obtained by the traditional method, given 



present computing capacity. A more promising direction 
is to use the transfer-matrix DMRG (TDMRG) method 
to handle a discretized imaginary time direction of 
finite size 25]. This should give good results at high 
temperature, and permits access to thermodynamic 
quantities (but not long-distance correlation functions). 
However, one would then have to abandon the naive 
chemical potential term used in the present work and 
add an imaginary, constant vector potential as is usual 
in Euclidean lattice studies [l|. 



D. Measurements 

For expectation values of observables defined on a sin- 
gle lattice site it is desirable to take measurements at the 
center of the lattice, using the single inserted site. Then 
for an observable O (for example jo) defined on the single 
site, we have 



(0) = 



E 

ijj'k 



4TP(x)h>4% = K po- 



ps) 



We calculate the correlation length £ in the system 
by looking at the decay of the two-point function. This 
serves as a useful cross-check to verify that the onset of 
a nonzero particle density occurs, as one would expect, 
close to fi — m at weak coupling. Unlike in Monte Carlo 
simulations, we do not have access to anything other than 
equal time correlations. 

To calculate the two-point function 

c ab {x - y) = (Mx)Mv)) - (M*)> (Mv)) > (24) 

we must store the operators for 4> a (x) and 4>b{y) until we 
are ready to calculate the two-point function 

(M x )Mv)) = E w%*lM x )]ii'[My)]kvWj k > 

■ii' jkk' 



and disconnected pieces 

(M*)) = T,^l*lM*)]<i>1>®k 



(25) 



(26) 



ii' jk 



(My)) = J2ri?kMy)hvri-l- (27) 

ijkk' 

So that the sum is over a complete set of states, it is 
simplest from a practical point of view if the two oper- 
ators come from different DMRG blocks [18(. Since our 
main interest in the two-point function is to determine 
the renormalized mass from the long-distance behavior, 
there is no reason to choose operators on the same block. 
Therefore we sacrifice measurements of the short-distance 
behavior to give better results at longer distances. 
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III. RESULTS AND COMPUTATIONAL 
CONSIDERATIONS 

The renormalized block Hamiltonians H are Hermitian 
matrices of size M 2 NxM 2 N. The main limitation on the 
accuracy and quality of the calculations presented here is 
the memory required to store these matrices (even when 
represented in a sparse format, this is over a gigabyte), 
and to a much lesser extent the additional computation 
time required by /i ^ 0. One "sweep" to the right or the 
left in the algorithm takes several hours of computer time 
for the largest truncations used here, although this can 
be less depending on the value of for fi = the Lanc- 
zos algorithm is used and performance is much better. 
It also depends on the number of eigenvectors required, 
making T^O prohibitively expensive given current com- 
putational resources. 

Computations were carried out with exhaustive 
searches of the parameter space. Numerical diagonaliza- 
tion routines from NAG were used; the superblock was 
stored as a sparse matrix. For fi = 0, the Hamiltonian is 
a real, symmetric matrix and the Lanczos method can be 
used. The current work, however, is primarily concerned 
with fi 0, for which it is easiest to use an implicitly 
restarted Arnoldi algorithm. For all the results shown, 
we tolerate a fractional error of at most 10 -8 in the nu- 
merically obtained Ritz vectors. 

The ground state energy is obtained directly when di- 
agonalizing the Hamiltonian. Other measurements are 
made as described in Sec. Ill Dl We wait until results 
change by no more than 1 part in 10 6 between sweeps 
before considering a quantity adequately converged. 

One naturally expects the charge operator Q to be self- 
adjoint so that in adding a term — u.Q to form an effective 
Hamiltonian we must have \x real. Unlike in calculations 
where the chemical potential is analytically continued, we 
have not made any changes to the original model. With 
enough memory and computer time we could work with 
H ^> m 2 , M and N very large and study nonzero particle 
densities far from the onset of condensation. There is no 
sign problem, merely the problem of representing a con- 
densate using the DMRG. However, it is accepted that 
the picture given by a finite truncation serves to capture 
the physics Q. 
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FIG. 1. Truncating the tower of bosonic states, with L — 40, 
m 2 = 1, A = 0.1 and M = N. The ground state energy is 
shown as a function of jj, for various consistently truncated 
bases. 



tower that must be truncated. One might choose to 
fix either the maximum number of bosons of any type 
{\h j) \i < n i j < n}: or the maximum number of bosons 
of either type {\i, j) < n}. It seems better to adopt 
the latter organization, as this is the truncation that bet- 
ter respects the global symmetry of the theory. Consider 
a field transformation 

4>i 0icos(0) + 2 sin(0); 
<p 2 -> <£icos(0) - 2 sin(0). 

To quantize this equivalent field theory, we can define 
new creation and annihilation operators in terms of the 
old ones that will obey the same canonical commutation 
relations 

a\(x) —> a\(x) cos(0) + a\(x) sin(0); 
a\(x) — > a\(x) cos(0) — a\(x) sin(0). 



A. Convergence with truncation size and 
arrangement 

Noting that storage considerations are a limiting fac- 
tor, we must consider the most reasonable approximate 
basis for the lattice sites. For simplicity, let us take the 
number of sites in the renormalized blocks M to equal 
those of the single site N. Finite-truncation effects are 
more noticeable when varying M, the renormalized block 
size. 

We choose a consistent way of truncating the bosonic 
states at each site as, in principle, there is an infinite 



Therefore if we use the truncated basis {\i, j) \ i+ j < n}, 
then we will always be able to reexpress the rotated state 
using the truncated basis whereas starting with the trun- 
cated basis j) \i < n,j < n} risks leaving us outside 
the span. Given this choice, Figure [T] shows how varying 
the number of states in the renormalized blocks affects 
the ground state energy. 

It is clear that, for a bosonic system, the DMRG con- 
verges quickly for zero particle density. It can also be 
reliably used to determine the onset of nonzero particle 
density. However, for large particle densities the trunca- 
tion is overwhelmed and accuracy is greatly diminished. 
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FIG. 2. Plot of the charge density jo as a function of \x at 
various L, for m 2 = 1, A = 0.1, and M = N = 28. There 
is no discernible finite-size effect for the spatial lattice size. 
It should be emphasized that for fi smaller than that shown 
here the charge density is exactly zero down to /i = 0. 



FIG. 3. Plot showing the ground state energy per site E/L 
as a function of /j, for m 2 = 1, A = 0.1, and M = N — 28 at 
various L. The finite-size effect is much more prominent here 
than in Fig. d 



B. Condensation and finite-size effects 



For A = 0.1 and various lattice sizes L, the particle 
density jo measured at the center of the lattice is shown 
in Fig. [2] These show a transition to a nonzero particle 
density as [x is increased, with no discernible dependence 
on L. On the other hand, the ground state energy per 
site - shown in Fig. [3]- has a very clear dependence on 
volume. The volume dependence is qualitatively similar 
for all values of the chemical potential studied. 

In Ref. [ll| it was noticed that more sweeps were re- 
quired for convergence close to the critical coupling. For 
fi close to condensation, a similar effect is observed. 



A=0.1 
X=0.2 
A=0.5 



C. Two-point function and phenomenological 
accuracy 

As A is increased from zero, the chemical potential 
at which the particle density becomes nonzero moves 
away from the free- field result. Figure 0] shows this for 
L = 40. We would expect that the condensation occurs 
when jj, > m^j, the renormalized mass, at weak coupling. 
The behavior of the two-point function is shown in Figure 
[5] for various A. We obtain the scalar masses nonpertur- 
batively from these results by a fit to the long-distance 
behavior C(x) = Ae~ mx j \fx. As is familiar from studies 
in lattice Monte Carlo, the "plateau" in the fitting errors 
is taken, to eliminate undesirable short-distance behav- 
ior (typically, the plateau is found at x s» 6). Finite- 
size effects were determined to be negligible by compar- 
ison with similar measurements taken with L — 10 and 
L = 20. Similarly, there was no discernible difference 



FIG. 4. The charge density jo is plotted as a function of /i, 
for m 2 = 1, L = 40, and M = N = 28 at various A. The 
phenomenologically expected zero charge density is observed 
for fi < fj,c- 

with results for a smaller truncation on the same lattice 
volume (L — 40 and N = M = 21). The measurements 
from these calculations agree within errors with the re- 
sults summarized in Table [U 

Taken together, these plots demonstrate that the 
present technique sidesteps the sign problem, and only 
gives a nonzero particle density in the phenomenologi- 
cally expected region /i > m^. As shown in Table HI we 
can separately estimate the value of the chemical poten- 
tial at which we see nonzero particle density from the 
data shown in Figure 21 and the renormalized mass from 
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FIG. 5. Example measurements of the two-point correlator 



C(x) for various A, m 
tion M = N = 28. 



1 with volume L = 40 and trunca- 



fits to data comparable to that of Fig. [SJ Taking the 
onset of nonzero particle density as occurring halfway 
between the last parameter choice where (jo) = and 
the first for which (jo) 7^ 0, the numbers are consistent 
for small A. 



IV. CONCLUSIONS 

We have shown that it is possible to use the DMRG 
to study bosonic, relativistic quantum field theories with 
two components in (1 + 1) dimensions at zero tempera- 
ture but nonzero chemical potential. The method cor- 
rectly captures the formation of a nonzero particle den- 
sity for large chemical potentials. It is our hope that it 
could serve as a useful numerical bridge from work in 
condensed-matter systems using microscopic Hamiltoni- 
ans to studies of relativistic field theories at finite density. 

Our discussion has considered various sources of er- 
ror, and how best to capture the nonperturbative physics 
with a finite truncation in boson number states. With 
improvements in numerical algorithms, and by tak- 
ing advantage of developments originally employed in 
condensed- matter contexts, the DMRG shows promise as 



an alternative to other numerical methods in relativistic 
field theory. With better computing resources, the re- 
sults described here can be extended to finite tempera- 
ture through the techniques described in Sec. Ill CI 

The applicability of the DMRG to bosonic systems will 
always be limited by the truncations and approximations 
involved: to a finite basis per site, to a finite volume, 
and to machine precision or worse in convergence of the 
Arnoldi algorithm. Fermionic systems do not suffer from 
the first of these three issues. Hence, as a step towards 
the long-term goal of studying QCD, one could revisit 
the massive Schwinger or Thirring models with nonzero 
chemical potential and at finite temperature [26| . 

Another innate limitation of DMRG is to one spatial 
dimension. Fortunately, much progress has been made 
by approaching the problem from a different direction. 
By identifying suitable tensor network states and corre- 
sponding renormalization procedures, systems in higher 
spatial dimensions can be studied. Such approaches in- 
clude projected entangled pair states [27J and the multi- 
scale entanglement renormalization ansatz [2^,[2^]. With 
these methods, the major issue to overcome is the large 
bond dimension needed to treat the system discussed in 
the present work; this will make adequate system sizes 
computationally too expensive for the time being [30| . 

The work in this paper has been carried out in the 
grand canonical ensemble. We can work in the canoni- 
cal ensemble by finding low-lying eigenvectors of the to- 
tal charge Q. This compels us to simultaneously diag- 
onalize H and Q. If we have a potential with interac- 
tions that allow a stable Q-ball to form (3j HH , then we 
would anticipate that such an inhomogeneity would ap- 
pear. This would allow the nonperturbative, relativistic 
study of nontopological solitons, even if the possibilities 
are limited by the lack of dynamics and the restrictive 
dimensionality. 

Perhaps one can also go the other way, from particle 
physics to condensed matter. A single-component Bosc- 
Hubbard model can be used to study nontopological soli- 
tons in nonrelativistic field theory; these may be thought 
of as Q-ball analogs [3J|. With such a simple model, it 
is possible to explicitly work in a sector of fixed charge 
by choice of basis [l6| , and then one simply diagonalizes 
the Hamiltonian as usual to find the ground state and 
excitations of the nontopological soliton. The interac- 
tions necessary to create such an object have not, to our 
knowledge, been studied with DMRG in either condensed 
matter or particle physics. 



A Mr /icon 

0.1 1.071 ±0.001 1.065 ±0.005 
0.2 1.114 ±0.001 1.123 ±0.003 
0.5 1.231 ±0.001 1.275 ±0.005 



TABLE I. Comparing renormalized masses and values of the 
chemical potential at which condensation occurs. For tur, the 
error quoted is the estimated error in the plateaued fit whereas 
for /i con it is half the spacing between parameter choices. 
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Appendix A: Matrix elements 

Here we give the matrix elements for a position-space 
discrctized two-component scalar field with an 0(2) sym- 



metry, an extension of the work in Ref. [ll|. We begin 
with the operators for the individual fields and their con- 
jugate momenta. On a given lattice site with basis states 
|m, n) it is straightforward to see, given ([S][7]) and requir- 
ing normalization, 



(m, n\ <f>\ |m', n') 
(m,n\(p2 \m',n') 
(to, n\ 7Ti |m', n') 
(to, n\ 7T2 |m', n') 



1 

71 
i 

71 

i 

71 

i 

71 



V TO — 


l^m— l,m ; 


+ Vto' - l<5 m , m 


-l 




\Jn — 


lo n -l,n' "1" 


- yV - 15„,„'-i 


5„ 




vto — 


l^m — 


- Vto' - l<5 m , TO 


-l 


(5 rii „/ and 



\/n- 15 n _1 „/ - V"' - l^r, 



(Al) 
(A2) 
(A3) 
(A4) 



(A5) 
(A6) 
(A7) 



To construct higher moments of the fields, we must write explicitly 

(to, n\ TT a n a \m! , n') = (m, n| 7r 2 + 7T 2 |m', n') , 
(to, n| <fi a 4>a \fn' i n') — (m, n\ cf)\ + </>?, |m', n') and 
(m, n| (0 a 0a) 2 |m', n') = (m, n| ^ + 2<^f^>| + ^ I™', n') . 

For the momenta components we then have 

(m,n\ tt\ \m',n') = - -\J(m - l)(m - 2)5 m -\, m '+\ + (2m - l)5 m>m / - \J (to' - l)(m' - 2)£ m+ i, m /_i <!>„,„' (A8) 
and 



(TO,n|7T2 |m',n') = - -\/(™ - - 2)<y n -i,n'+i + (2n - 1)<5„,„' - yj (n' - l)(n' - 2)<5„ + i, n /-i 

Finally, for the field components we have 
1 



(A9) 



(m, n\ (f>1 \mf, n') 
(m,n\ 4>l \m',n) 
(to, n\ 4>\ \m', m!) 



\J{m- l)(m - 2)<y m _i, m /+i + (2m - l)S m , m > + V Jjn' ~ l)(m' - 2)<J m+ i, m <_i <5„ >n /, (A10) 



V>- l)(n-2)<J n _i, n /+i + (2n - + VV - - 2)<5 n+ i,„< 



l r 



\J(m- l)(m - 2)(to - 3) (to - 4)5 m _i, m '+3 



(All) 
(A12) 



l,m' + l 



+ 4to-\/ (to — l)(m — 2)5 m , m '+2 - 6 V (m - l)(m - 2)5 m _ 

- 6\/(m' - 1)(to' - 2)5 m +i, m '-\ + 3(2to 2 - 2m + l)£ m , m < + 4mV(™' - 1)(to' - 2)* 



m+2 



+ V 7 ~ !)( m ' ~ 2 )(to' ~ 3 )( TO ' ^ 4)(5 m+3 , m '_i 



(TO,n| <t4 \m',n') = - y/(n-l)(n- 2)(n - 3)(n - 4)<5 n _i, n / +3 



(A13) 



+ 4nv/ (n - l)(n - 2)<5 n ,„< + 2 - 6a/ (n - l)(ra - 2)<S n -i, n '+i 



- 6V(n'-l)(n / -2)* B+ i, n /_i + 3(2 ? i 2 - 2n + l)£„, n , + 4«'vV - !)(»' - 2)<Wa,n' 
+ \/V - l)(n' - 2)(n' - 3)(n' - 4)5„ +3 ,n'-l <5m,m'- 
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